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ABSTRACT 

We present a possible star formation and chemical evolutionary history for two early- 
type galaxies NGC 1407 and NGC 1400. They are the two brightest galaxies of the 
NGC 1407 (or Eridanus-A) group, one of the 60 groups studied as part of the Group 
Evolution Multi- wavelength Study (GEMS). 

Our analysis is based on new high signal-to-noise spatially resolved integrated 
spectra obtained at the ESO 3.6m telescope, out to ~ 0.6 (NGC 1407) and ~ 1.3 
(NGC 1400) effective radii. Using Lick/IDS indices we estimate luminosity-weighted 
ages, metallicities and a-element abundance ratios. Colour radial distributions from 
HST/ACS and Subaru Suprime-Cam multi-band wide- field imaging are compared 
to colours predicted from spectroscopically determinated ages and metallicities using 
single stellar population models. The galaxies formed over half of their mass in a single 
short-lived burst of star formation (^ 100 MQ/year) at redshift z ^ 5. This likely 
involved an outside-in mechanism with supernova-driven galactic winds, as suggested 
by the flatness of the a-element radial profiles and the strong negative metallicity 
gradients. Our results support the predictions of the revised version of the monolithic 
collapse model for galaxy formation and evolution. We speculate that, since formation 
the galaxies have evolved quiescently and that we are witnessing the first infall of 
NGC 1400 in the group. 

Key words: galaxies: individual: NGC 1407, NGC 1400 - galaxies: formation - 
galaxies: evolution - galaxies: stellar content - galaxies: abundances 



1 INTRODUCTION 

The physical formation and evolution processes taking place 
at various locations within a galaxy leave their fossil imprint 
in the stars. These chemodynamical imprints are the regu- 
lar features observed in the stellar populations, photome- 
try, internal kinematics and structure of early-type galaxies. 
Furthermore, these features are found to correlate defining 
scaling relations of remarkable tightness, such as the colour- 
magnitude relation and the fundamental plane. This class 
of galaxy, therefore, is one of the best guides to compare 
observations and theoretical model predictions. 

Theoretical scenarios that describe the formation and 
evolution of early-type galaxies include: a revised version of 
monolithic collapse, the dissipative (wet/gas-rich) and dissi- 
pationless (dry/gas-poor) merger alternatives of hierarchical 
clustering. 
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In the classical models 
lapse jEggen et al. 1963: 



of m onolithic 

Larson 1974al : 



col- 

■ _^ ■ iLarson 19751 : 

ICarlberg 19841 : lArimoto fc Yoshii 1987h . stars form in 
all regions during a rapid collapse and remain in their 
orbits, whereas the gas dissipates to the centre of the 
galaxy. The shape and ellipticity of stellar orbits are 
expected to depend on the initial angular momentum of 
the collapsing protogalactic gas cloud and on the effect 
of gas viscosity in redistributing angular momentum, 
towards outer radii, during the collapse. The sinking gas 
is continuously enriched by evolving stars. As a conse- 
quence, stars formed in the centre are predicted to be more 
metal-rich than those born in the outer galaxy regions. 
The efficiency of this process is proportional to the depth 
of the galaxy potential well, so that a strong correlation 
between metallicity gradient and galaxy mass is expected. 
The evolving stars are also responsible for the a-element 
enrichment of the interstellar medium (ISM). To reproduce 
the supersolar abundance of a-elements with respect to 
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the iron-peak elements, observed in early-type galaxies, 
the models require that the collapse and star forming 
process occurred on timescales < 1 Gyr, which will produce 
null or very small age gradien ts (|Arimoto fc Yoshii 19871 : 
iMatteucci 1993 : iThomas 19991 ). Recently, feedback pro- 
cesses such as supernova-driven galactic winds have been 
re-considered in more detailed numerical models of dis- 
sipative collapse ( Larson 1974bl: lArimoto fc Yoshii 1987 ; 



Gibson 19971: iMartinelli et al. 1998 : IChiosi fc Carraro 2002 
Kawata fc Gibson 20031 : iPipino et al. 200e 

Pipino et al. 200*71 )! It is suggested that they may cre- 



ate and shape the abundance gradients. Galactic winds 
initiate when the energy injected by supernovae explosion 
balances the gravitational binding energy of the ISM. The 
winds evacuate the gas, eliminating the essential fuel for 
any future star formation. The shallower local potential 
well of the external galactic regions supports the early 
development of winds with respect to the central regions. 
Therefore, in the central regions the star formation and the 
chemical enrichment last longer. Dissipation of gas towards 
the galaxy centre and a time-delay in the occurrence of 
galactic winds cooperate in steepening any metallicity 
gradient. In summary, dissipative collapse models predict: 
i) very steep negative metallicity gradients, that strongly 
correlate with galax y mass (jChiosi fc Carraro 20021 : 
iKawata fc Gibson 20031 ). ii) small or null age gradients, as a 
consequence of the short timescales involved in the collapse 
and star formation processes, iii) a-element enhancement 
yadients that can b e either positive, n egative or null 
Martinelh et al. 19981 : iPipino et al. 20061 ). due to the 
"inter-play between local differe nces in the star formation 
timescale and gas flows" ((Pipino et al. 20071 ). iv) isophotal 
shape and ellipticity proportional to the protogalactic 
initial angular momentum and its radial redistribution, v) 
no peculiar internal kinematic structure. The difference 
between the original and revised dissipative collapse models 
focuses on the assembly of the initial gaseous material. The 
revised view accounts for the possibility of either a unique 
primordial cloud or the coalescence of many gaseous clumps 
without any preexisting stars. Otherwise, the physical 
processes in galaxy formation are similar to that described 
above. 

In the hier archical clustering sc e nario of galaxy 



formation (eg. 


Kauffmann et al. 19931: Cole et al. 1994; 


BauKh et al. 1996 


: iKauffmann fc Chariot 1998 ; 


De Lucia et al. 20061). elliptical galaxies are produced 



by merging events. Early numerical simulations of galaxy 
mergers predicted contradictory r adial variatio ns of stellar 
properties in merger remnants. IWhite (1980l ) suggested 
that dissipationless mergers cause a flat tening of metallicity 
gradients, whereas Ivan Albada (1983 ) argued that the 
position of the stars in the local potential are preserved 
by violent relaxation and the gradients in the progeni- 
tors are only affected mo derately. Further simulations of 
iBarnes fc Hernquist (199ll ), considered the hydrodynamical 
physics of the gas and found that during a dissipative 
merger a significant fraction of the progenitors' gas tends to 
migrate toward the centr al regions of the merger remnant. 
iMihos fc Hernquist (1994 ) showed that the gas accumulated 
in the central regions may trigger a local secondary burst 
of star formation. Galaxy merger-induced star formation is 
predicted to produce metallicity and age gradients. Positive 



or negative a-element enhancement gradients may also 
originate, depending on the original enhancement of the 
gas and the duration of t he burst (Thomas et al . 19991 : 
iThomas fc Kauffmann 19991 ). iBekki fc Shiova (1999i ) sim- 
ulated dissipative mergers of two gas-rich disc galaxies of 
varying mass to form an elliptical galaxy remnant. Their 
models show that metallicity gradients in merger remnants 
are shallower with respect to those predicted by dissipative 
collapse model s and only weakly dependent by the remnant 
galaxy mass. iKobavashi (2004 ) simulated the formation 
and chemodynamical evolution of elliptical galaxies in a 
cold dark matter cosmology, focusing on internal metallicity 
gradients. The models ranged from a monolithic collapse 
scenario, seen as the assembly of tens of gas-rich subunits 
at high redshift, to dissipative and dissipationless mergers 
of equal-mass galaxies at low redshift (which she denoted 
"major mergers"). The result was that galaxies of a 
given mass had steep metallicity gradients if formed by 
collapse and shallower gradients if formed by major merger 
(dissipative or dissipationless). 

Recently, combining t he Millennium co smolog- 
ical N-body simulation llSpringel et al. 20051 ) with 
semi-analytic models, Lucia et al. (20061 ) and 

Lucia fc Blaizot (20071 ) focused on the distinction 
between the formation and assembly time (redshift) of 
elliptical galaxies and brightest cluster galaxies (BCGs) in 
the hierarchical clustering scenario. The models make use 
of the prescription of feedback from a central active gala ctic 
nucleus (AGN) and supernova (see ICroton et al. 200 ^ to 
prevent excessive galaxy-mass growth from cooling flows. 
Their model predicts a mass-dependent evolutionary history 
("down sizing" scenario), with more massive galaxies (e.g 
BCGs, BGGs; brightest group galaxies) forming a large 
fraction of their stars (50 — 80 per cent) in progenitor 
systems at redshift ~ 2.5. The progenitors then assemble 
into a single final object at z ~ 0.8. The formation epoch 
of less massive galaxies is predicted to be z ~ 1.9 and 
the assembly time at redshift ~ 1.5. The peak of the star 
formation history of massive galaxies is at redshift ~ 5, 
and it progressively moves towards lower redshifts for less 
massive systems. Similarly, the duration of the star forming 
episode is predicted to last longer in low mass galaxies. 
Mergers are thought to induce bursts of star formation. 
Therefore, the star formation histo ry of a galaxy present s 
the "bursty" behaviour described in lDe Lucia et al. (20061 ). 

The kinematics and internal structure of a merger 
remnant are a consequence of the progenitors' mass ra- 
tios and the ef fi ciency of the dissipat ive process (e.g. 
iNaab et al. 20061 ). iNaalj fc Burkert (20031 ) show that a dis- 
sipationless binary merger of equal-mass disk galaxies leads 
to a slowly rotating and anisotropic supported system 
{(v/a)* < 0.4) with preferentially boxy isophotes and signif- 
icant minor-axis rotation. Whereas, the merger of unequal- 
mass disk galaxies produces a rotationally supported ellip- 
tical galaxy with a small amount of minor-axis rotation 
{{v/a)* — 1.2; suggesting isotropic oblate rotators) and 
disky isophotes. In dissipative mergers the gas is found to 
settle at the galaxy centre deepening the potential well and 
making it more axisymmetric. The isophotal shape of an 
equal-mass merger is then strongly affected by the gas, since 
most of the stars are gravitationally induced to move from 
boxy to more elliptical or even disky orbits. Similarly, the 
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disky isophotal shape of an unequal-ma ss merger remnan t 
is weakly affected by the presence of gas (|Naab et al. 200a ) . 
Furthermore, in simulat ions of an equal-m ass merger with 
induced star formation (|Bekki fc Shiova~1997. ) gradual gas 
dissipation forms an elliptical galaxy with disky isophotes, 
whereas a merger with rapid star formation produces both 
boxy and disky isophotes depending on the viewing angle 
of the observer. The merger of unequal-mass gas-rich disk 
galaxies is found to exhaust a large amount of the gas ow- 
ing to moderate sta r formation and ending in the forma- 
tion of an SO galaxy (|Bekki 19981 ) . In summary, hierarchical 
galaxy formation models predict: i) shallow negative metal- 
licity gradients, with little host mass dependence, ii) age gra- 
dients and iii) a-element enhancement gradients, that can 
be either positive or negative, due to the duration and loca- 
tion of the merger-induced star formation and the original 
abundance pattern of the gas, iv) isophotal shape and el- 
lipticity depending on the progenitors mass ratio and the 
efficiency of gas dissipation, v) peculiar internal kinematic 
structures. 

The contrasting predictions of the galaxy formation 
models underline the lack of complete understanding of 
the physical processes behind the formation of galaxies. 
Observationally, we now have the ability to obtain spatially 
resolved high signal-to-noise integrated spectra out to large 
galactic radii (and hence larger mass fractions). This allows 
us to consider the physical mechanisms acting locally and 
how their properties vary with radius (i.e. investigating be- 
yond the central regions), proving invaluable constraints on 
the processes of galaxy formation and evolution. Recently, 
an increasing number of works have focused their attention 
on stellar population radial profiles of galaxies in dif- 
ferent environmen ts 
Mehlert et al. 200 j; 



Name 


Type 


PA 


re 


Ms 








(degrees) 


(arcsec) 


(mag) 


(mag) 


NGC 1407 


EO 


35 


72" 


-21.22 


-24.77 


NGC 1400 


SO 


40 


27" 


-19.95 


-23.74 



(e.g. Kobavashi fc Arimoto 199S 



ISanchez-Blazquez et al. 200(: 

20071 : iBroueh et al. 2007 



Sanchcz-Blazguc z et al 
Reda et al. 20071) . The results underline a common 
behaviour for the stellar population radial profiles of the 
studied galaxies: early- type galaxies are often characterised 
by strong metallicity gradients, shallow age gradients and 
null or statistically insignificant a-element enhancement 
gradients. 

Following from ISpol aor et al. (2008al ) (hereafter Pa- 
per I) we probe a possible star formation and chemical 
evolutionary history of the early-type galaxies NGC 1407 
and NGC 1400. W e accomplish the task using the spectra l 
Lick/IDS indices (|Faber et al. 19851 : IWorthev et al. 19941 ) 
to determine the radial profiles in age, metallicity, [Z/H], 
and a-element abundance ratios, [a/Fe], out to ~ 0.6 
(NGC 1407) and ~ 1.30 (NGC 1400) times the galaxies' ef- 
fective radii (re; the radius within which half the galaxy light 
is contained). Central values and gradients in age, [Z/H] 
and [a/Fe] are estimated and used as proxies of the quan- 
tity, velocity and duration of gas dissipation, star formation 
and possible merger events during galaxy formation. Fur- 
thermore, we have considered the results from the spatially 
resolved radial kinematics and surface photometry analysis 
performed in Paper I. We compared colour radial profiles 
with colours predicted from spectroscopically derived ages 
and metallicities using single stellar population models. 

This paper is organised as follows. In Section 2 we de- 
scribe the two sample galaxies. In Sections 3 and 4 we de- 
scribe the observations and the relevant data reduction. Sec- 



Table 1. Galaxies properties. Type: morphological type extracted 
from NASA/IPAC Extragalactic Database (NED); PA: position 
angle of the major axis; r^: effective radius from HyperLEDA; 
Ms, M/f: absolute B-band and K-band magnitudes obtained 
from HyperLEDA and 2MASS, respectively. 



tion 4 presents the spatially resolved stellar population pa- 
rameters. In Section 5 observed and predicted colour index 
radial profiles are examined. In Section 6 we summarise our 
results, and discuss the possible star formation and evolu- 
tionary histories of NGC 1407 and NGC 1400. In Section 7 
conclusions are presented. 



2 THE DATA SAMPLE 

The NGC 1407, or Eridanus-A, group is one of the 
60 groups studied as part of the Group Evolution 
Multi-wavelength Study ('GEMS: lo"smond fc Ponman 20041 : 
iForbes et al. 20061 ). The study was designed to probe the 
evolution of groups and their member galaxies. NGC 1407 
is the brightest group galaxy of the NGC 1407 (or Eridanus- 
A) group. NGC 1400 is the second brightest group galaxy 
and it is characterised by a peculiar low velocity. A com- 
plete description of the data sample is presented in Paper I. 
Table [1] reports some of the galaxies' properties. 



3 OBSERVATIONS AND DATA REDUCTION 

Spectral observations were performed with the ESO Faint 
Object Spectrograph and Camera (EFOSC2) mounted on 
the ESO 3.6m telescope at La Silla Observatory, Chile. Ta- 
ble [2] summarises the instrumental configuration adopted 
during the observing run. Lick/IDS and spectrophotometric 
standard stars were taken at the parallactic angle, the lat- 
ter with a 5" wide slit. The Lick/IDS stars are also used as 
velocity standards. 

Data reductions were carried out with IRAfQ using the 
procedures adopted in Paper I. In Fig. [1] is presented the 
spectrum of NGC 1407 and NGC 1400, obtained from co- 
adding the individual ID spectra. 

We spatially resolved 68 and 82 apertures along the 
observing axes for NGC 1407 and NGC 1400, respectively. 
The spatial width (i.e. the number of CCD rows binned) for 
each extracted aperture increases with radius to achieve a 
signal-to- noise ratio of ^ 30 at 5000 A. We reached a 
radial extent of ~ 0.56re (~ 4.11 kpc) for NGC 1407 and 
~ l.SOre (~ 3.58 kpc) for NGC 1400. 



^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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F.SO-'^ fim 


Spectrograph 


EFOSC2 


Chip size (pixels) 


2048 X 2048 


Grism resolution FWHM 


7.8 A 


Pixel scale 


0.314" 


Binning 


2x2 


Slit width 


1.2" 


Wavelength coverage 


4320 - 6360 A 


Seeing FWHM 


~ 1" 


Exposure time 


3 X 1200 s 


NGC 1407 Slit P.A. 


44° 


NGC 1400 Slit P.A. 


42° 



Table 2. Spectral observing parameters. 
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Figure 1. The spectrum of NGC 1407 and NGC 1400, obtained 
from co-adding the individual ID spectra. The flux scale is ar- 
bitrary. The shadowed regions indicate the central band of the 
Lick/IDS indices used in the stellar populations analysis. 



4 SPATIALLY RESOLVED STELLAR 
POPULATIONS 

The wavelength range of 4300-6300A covered by our obser- 
vations includes 16 Lick/IDS indices jWorthev et al. 1994 ). 
These include one Balmer index (H/3), three Magne- 
sium indices (Mgi, Mg2 and Mgi,), Calcium (Ca4455), 
Sodium (NaD), Titanium Oxide (TiOi) indices plus nine 
Iron indices Fe4383, Fe4531, Fe4668 (referred as C4668), 
Fe5015, Fe5270, Fe5335, Fe5406, Fe5709 and Fe5782. The 
indices were measured using the method described in 
iProctor fc Sansom (20021 ). iTrager et al. (199 81) index defini- 
tions were adopted in the analysis. The calibration proce- 
dure of our indices measurements to the Lick/IDS system is 
described in Appendix 1X1 

Fully corrected Lick/IDS line-strength indices, ex- 
tracted from the aperture spectra along the observing axis of 
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Figure 2. The average deviation in units of error (i.e. x) for 
the 16 observed Lick/IDS indices. Indices are listed in increasing 
wavelength order. Error bars represent the rms scatter in the 
deviations. Deviations after the removal of the poorly fitted Mgi 
and NaD indices (open circles) and of the indices omitted by the 
iterated 3cr-clipping routine (see text for details). 



both the galaxies, are available electronically. In Figures [B II 
and IB2I the Lick/IDS index values are plotted against the 
galactocentric radius of NGC 1407 and NGC 1400 respec- 
tively. Logarithmic gradients of the indices were measured 
by applying the linear regression analysis described in Ap- 
pendixE The resuhs are reported in Tables[Bl](NGC 1407) 
and m (NGC 1400). 



4.1 Stellar population model fitting procedure 

The observed Lick/IDS indices have been used to de- 
rive luminosity- weighted log(Age), total metallicity [Z/H] 
and a-abundance ratio [a/Fe]. The a-abundance ratio is 
parametrised by [E/Fe]; the latter quantifies the enhance- 
ment from the a-elements (N, O, Mg, Na, Si, Ti) with re- 
spe ct to the Fe-peak ele ments (Cr, Mn, Fe, Co, Ni, Cu, Zn; 
see iThomas et al. 20031 for de tails) . 

Briefly, the technique jProctor fc Sansom 20021 ) in- 
volves the simultaneous comparison of as many observed 
indices as possible to single stellar population (SSP) mod- 
els. The best fit is achieved throughout a statistical x^- 
minimisation technique, which reduce the deviations be- 
tween the observed and the modelled values as a fraction 
of the index errors, i.e. . The strength of the method is 
that it works with as many indices as possible in order to 
break the age- metallicity degeneracy that affects each index 
differently . 

In this work, we use I Thomas et al. (20031 ) SSP models. 
Models are provided with: [a/Fe] = —0.3, 0.0, 0.3, 0.5; ages 
from 0.1 to 15 Gyr (-1 log(age) 1.175, in steps of 0.025 
dex); [Z/H] from -2.25 to 0.8 in steps of 0.025 dex. 

We began the model fitting process obtaining the best 
fit to all 16 observed indices. The residuals are expressed in 
terms of index errors y. We found that Mg i and NaD indices 
are poorly fit by the lThomas et al. (20031 ) SSP models. In- 
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deed, it is well known that interstellar absorption aflects the 
strength of the NaD index. The Mgi index shows a small rms 
deviation but has an average deviation value greater then 
3a. This deviation is not caused by the 5102A sky emission 
line, however a large average deviation could be associated 
with poor flux calibration. In fact, the Mgi and Mg2 indices 
are known to be sensitive to flux calibration corrections. 
We detected some residuals from the 5577A sky emission 
line, but they did not affect any nearby Lick/IDS absorp- 
tion features. We performed a new fit removing the Mgi and 
NaD indices and applying an iterated 3 cr-clipping function 
l|Proctor et al. 20041 : iProctor et al. 20051 ) for the remammg 
indices. This lead to a better fit with improved residuals 
(Fig. 2). 



4.2 Central stellar population parameters 

The central stellar population parameters of NGC 1407 and 
NGC 1400 were measured by averaging the values of all the 
apertures within a radius of re/8. Central values of age, total 
metallicity, [E/Fe] ratio and velocity dispersion are reported 
in Table |3] together with values from previous studies. 

We find that our measurements of NGC 1407 
are in good agreement with literature values wit h the 
exception of the results of iDenicolo et al. (20051 ) and 
omas et al. (20051 ) ■ The former found a younger cen- 
tral age (2.5 Gyr), a more metal rich stellar popula- 
tion ([Z/H]=0.67) and less enhancement from a-elements 
([E/Fe]=0.18). The reason for this disagreement may be as- 
sociated with the em ission line corrections they applied. In 
iThomas et al. (20051 ) the central age is younger (7.4 Gyr) 
than ours, whereas the [Z/H]=0.38 and [E/Fe]=0.30 values 
are in good agree ment. For NG C 1400 our values are in goo d 
agreement with iHowell (20051) and Irhomas et al. (2005h . 
but only partiallv with lSarr et al. (2007h . In the latter, the 
central age is younger, the [E/Fe] ratio is in good agreement, 
but the total metallicity value is higher than our value. Thus, 
we find good agreement with the majority of literature val- 
ues. 

In summary, we find that the central regions of both 
the galaxies present old 10 Gyr), metal-rich stellar pop- 
ulations with supersolar abundances of a-elements. 



4.3 Stellar population radial profiles 

Stellar population parameters extracted from the spatially 
resolved apertures along the observing axes of NGC 1407 
and NGC 1400 have been plotted in Fig. O The logarithmic 
gradients of the radial profiles are measured using the same 
method adopted for Lick/IDS line-strength radial profiles 
(see Appendix [B} , beyond the seeing limit of 1 arcsec; the 
results are summarised in Table U The quoted errors are the 
standard errors on the slope. 

The radial age profile of NGC 1407 indicates a con- 
stant old age for the spatially resolved regions of the galaxy. 
The measured age gradient is not statistically significant: 
—0.00 ± 0.02 dex per radius dex. The region inside 5 arc- 
sec was analysed to search for variations in the age profile 
likely to be associated with the possible kinematically decou- 
pled core (Paper I). No significant substructure was found 
that could not be due to the intrinsic scatter of the data 



points. The measured age gradient of NGC 1400 is also not 
statistically significant: —0.02 ± 0.02 dex per radius dex. 
Previous works on age gradients of early-type galaxies in 
different environments (see Table |4]) are consistent with 
zero-slope or very shallow gradients. iMehlert et al. (20031 ) 
found no significant age gradients for 91 per cent of their 
sample of 35 early-type gal axies in the Coma cluster. 
Sanchez-Blaz guez et al. (20061 ). in a study of 76 early- type 
galaxies situated in low- and high-density environments 
did not find any overall significant age gradient in the 
sample. iReda et al. (20070 also found zero-slope age gradi- 
ents for 12 isolated early- type galaxies. In a recent study, 
ISanchez-Blazguez et al. (20071 ) found a shallow mean age 
gradient for 11 early- type galaxies covering a wide range 
in luminosity and environment. 

Both NGC 1407 and NGC 1400 galaxies present 
metallicity radial profiles similar in shape and 
slope. We measured gradients of —0.38 ± 0.04 
and -0.47 ± 0.04 dex p er radius dex, respectively. 
iKobavashi fc Arimoto (19991 ) derived a mean metallicity 
gradient from absorption-line strengths literature values 
of 80 earl y-type galaxies of —0.30 ± 0.15 dex per ra- 
dius dex. ISanchez-Blazguez et al. (20061 ) found a metal- 
licity gradient of —0.21 ± 0.02 dex per radius dex 
for low-density environment elliptical galaxies and 
—0.33 ± 0.06 dex per radius dex for high-density en- 
vironment elliptical galaxie s. Similar negative metallicity 
values were found in ISanchez-Blazguez et al. (2007 ) 
(-0.33 ± 0.07 dex per radius dex) and ' Reda et al. (20071 ) 
(—0.25 ± 0.05 dex per radius dex) for a sample of isolated 
galaxies. Our literature search (Table |4]) finds a range of 
variation for the measured metallicity gradients: i.e going 
from a minimum of grad [Z/H] — -0.16 dex per radi us dex 
(in the Coma cluster galaxies; iMehlert et al. 2( )03') to a 
maxim um of —0.34 dex per radius dex (for BGGs and 
BCGs; iBrough et al. 20071 ). 

Scatter along the fit is seen in both the [E/Fe] pro- 
files, but no statistically significant gradient was mea- 
sured. Statistica lly no significant rad ial gradients of [E/Fe] 
were found by IMehlert et al. (20031 ) with an average of 
0.05 ± .05 dex per radius dex. Null mean gradients were 
foun d bvlSanchez-Blazguez et al. (20071 ) , iReda et al. ( 20071 ) 
and iBrough et al. (20o'7i r Given the absence of significant 
gradients in the a-abundance ratio values, those derived 
from the central regions (r ^ re/8; Table [3| are represen- 
tative of the whole galaxy. Of particular interest is that 
both galaxies are characterised by a significant supersolar 
abundance of a-elements, [E/Fe]o = -1-0.41 and -1-0.33 re- 
spectively, suggesting a short star formation time scale. The 
[E/Fe] radial profiles, as well as the metallicity radial profiles 
and age, did not present any evidence substructure related 
to the possible kinematically decoupled core in NGC 1407 
(Paper I). 

In summary, we find that the stellar population radial 
profiles of both NGC 1407 and NGC 1400 galaxy are sim- 
ilar in shape and slope, suggesting a similar formation and 
evolution scenario. The radial age and [E/Fe] profiles are 
flat, with no statistically significant gradients measured. The 
metallicity radial proflles present significant negative gradi- 
ents, with an average value steeper than previous literature 
works on metallicity gradients of early-type galaxies. 
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Ageo 


[Z/H]o 


[E/Fc]o 






(Gyr) 


(dex) 


(dcx) 


(km s-i) 


NGC 1407 


This work 


12.0 ± 1.1 


0.29 ± 0.08 


0.41 ± 0.05 


278.9 ± 4.3 


Humphrev & Buote C2006') 


12.0 ± 2.0 


0.35t ± 0.06 


0.33 ± 0.02 


259.7 ± 3.7 


Cenarro et al. (2007) 


11.9 ± 1.4 


0.38 ± 0.04 


0.33 ± 0.01 


293.5 ± 4.1 


Cantiello et al. C2005~) 


11.0 ± 1.0 








Howell ('2005) 


9.5 ± 2.2 


0.56 ± 0.07 


0.30 ± 0.04 


296.0 ± 4.0 


Zhane et al. (2007) 


9.2 ± 0.7 








Thomas et al. (2005) 


7.4 ± 1.8 


0.38 ± 0.04 


0.30 ± 0.02 


259.7 ± 3.7 


Denicolo et al. ('20051 


2.5 ± 0.6 


0.67 ± 0.05 


0.18 ± 0.05 


265.0 ± 17.0 


NGC 1400 


This work 


13.8 ± 1.1 


0.25 ± 0.06 


0.33 ± 0.04 


254.4 ± 3.1 


Howell ("2005) 


14.2 ± 3.0 


0.31 ± 0.13 


0.35 ± 0.04 


285.0 ± 5.0 


Thomas et al. (2005) 


15.0 ± 1.6 


0.23 ± 0.03 


0.34 ± 0.01 


256.2 ± 3.0 


Barr et al. f2007) 


9.57 ± 1.1 


0.66t ± 0.03 


0.34 ± 0.01 


215 ± 7.0 



Table 3. Central (r < re/8) stellar population parameters for NGC 1407 and NGC 1400 from this work and from previous studies. The 
errors are on the mean. The sym bol f indicates [Z/H ] values that were not provided in the original works, but that we recovered using 
the [Z/H] = [Fe/H] + 0.94[o/Fe] dTraeer et al. 2000l) . 
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Figure 3. Spatially resolved radial profiles for the NGC 1407 (left panel) and NGC 1400 (right panel) derived stellar population 
parameters: log(Age), total metallicity [Z/H] and Q-elements abundance ratio [E/Fe]. Filled red squares represent values from the 
spatially resolved apertures extracted along the observing axis in the negative binning direction; open circles are the same in the positive 
binning direction (see text for details). The approximate extent of the seeing disc is identified by a dashed line. The error bars express 
the median error in the derived parameters. Gradients (black thick lines) are measured excluding the points inside the region affected 
by the seeing. The quantities are measured in logarithmic radius scaled by the effective radius of the galaxy. 



5 OBSERVED AND PREDICTED COLOUR 
INDEX RADIAL PROFILES 

In Paper I we presented the results of a spatially re- 
solved surface photometry analysis of NGC 1407 and 
NGC 1400. The study was performed using HST/ACS and 
Subaru Suprime-Cam multi-band imaging data. The high 
instrumental resolution of HST/ACS allowed us to sample 
NGC 1407 over radii ~ 0.1 - 100 arcsec (~ 0.01 - 10.20 kpc) 



in B and I band. The wide-field of Subaru Suprime-Cam al- 
lowed us to sample NGC 1407 and NGC 1400 over radii 
-5-100 arcsec (~ 0.51 - 10.20 kpc) in g', r' and i' Sloan 
band. The lost of the inner ~ 5 arcsec is caused by saturation 
within the CCD (see Paper I for details). 

In Fig. |4] we super-imposed the colour predic - 
tions from the SSP models of [Thomas et al. (20031) 
( http://www-astro. physics. ox. ac.uk/~maraston/SSPn/colors/| 
to the observed colour profiles. The models provides colours 
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N 


Type & Environment 


grad Age 


grad [Z/H] 


grad [E/Fe] 








(dex per dex) 


(dex per dex) 


(dex per dex) 


NGC 1407 


1 


E; group 


-0.00 ± 0.02 


-0.38 ± 0.04 


0.02 ± 0.02 


NGC 1400 


1 


SO; group 


-0.02 ± 0.02 


-0.47 ± 0.04 


0.04 ± 0.03 


Averages 












This work 


2 


E, SO; group 


-0.01 ± 0.02 


-0.42 ± 0.04 


0.03 ± 0.02 


Gorgas et al. f 1990) 


15 


Es, SOs, BCGs; groups, clusters 




-0.22 ± 0.10 




Fisher et al. (1995) 


16 


BCGs, Es; groups, clusters 




-0.25 ± 0.10 




Kobavashi & Arimoto (1999) 


80 


Es; groups, clusters 




-0.30 ± 0.15 




Mehlert et al. C2003) 


35 


Es, SOs; Coma cluster 


0.43 ± 0.46 


-0.16 ± 0.12 


0.05 ± 0.05 


Sanchcz-BlazQuez et al. f2006) 


61 


Es, SOs; field, group, Virgo cluster 


0.08 ± 0.02 


-0.21 ± 0.02 




Sanchcz-Blazquoz ot al. (2006) 


15 


Es, SOs; Coma cluster, Abell clusters 


0.03 ± 0.07 


-0.33 ± 0.06 




Sanchcz-Blazqucz ct al. (2007) 


11 


Es, SOs; field, groups, Leo group, Virgo cluster 


0.16 ± 0.05 


-0.33 ± 0.07 


0.00 ± 0.05 


Brough et al. (2007) 


6 


BCGs, BCGs; groups, clusters 


0.01 ± 0.04 


-0.34 ± 0.08 


0.06 ± 0.02 


Reda et al. (2007) 


12 


Es, SOs; isolated 


0.04 ± 0.08 


-0.25 ± 0.05 


-0.03 ± 0.02 



Table 4. Stellar population radial gradients for NGC 1407 and NGC 1400. Average values of age, [Z/H] and [E/Fe] gradients in early- 
type galaxies obtained in other works. N: number of galaxies considered in the study. Type & Environment: morphological type and 
environment of the studied galaxies. Grad Age: age gradient, —^^j^^y-- Grad [Z/H]: total metallicity gradient, ^j^pj- Grad [E/Fe]: 
a-elements abundance ratio gradient, -^^^t^- The quoted errors are the standard errors on the slope. 



as a function of age and metallicity. Therefore, using our 
spectroscopically measured ages and metallicities (see 
Section 4.3) we created radial profiles of predicted colours. 

The predicted g'— i', g'— r' and r'— i' profiles of both the 
galaxies required ad hoc vertical offsets to match the ob- 
served colours: i.e. —0.04, —0.07 and -1-0.03, respectively. 
We compared observed and predicted colours by measur- 
ing the radial gradient of the profiles. The results are re- 
ported in Table [5] The predicted NGC 1407 B-I pro- 
file presents a steep gradient that is not seen in the ob- 
served colours; whereas small differences are found between 
the Sloan colour profile grad ients. In a similar analysis on 
NGC 821, iProctor et al. (2005 ) also found it necessary to 
apply offsets in order to match observations. These offsets 
might be explained with effects of internal extinction, al- 
though systematic effects (e.g. SSP modelling errors) may 
be present. 

To better investigate the difference in slope we calcu- 
lated the ratio AC/A[Z/H] (see Table E}. In particular, 
{AC / A[Z/ H])obs is the ratio between the observed colour 
index gradient and the gradient in metallicity that we found 
from the Lick/IDS indices analysis (see Table In the 
[AC / A[Z/ H])ssp ratio, AC expresses the colour index dif- 
ference estimated from the SSP model assuming a constant 
age (12 Gyr for NGC 1407 and 13.8 Gyr for NGC 1400; see 
Table[3l) and metallicity values of [Z/H] = 0, [Z/H] = -1-0.38 
for NGC 1407 and [Z/H] = 0, [Z/H] = +0A7 for NGC 1400; 
consequently, A[Z/H] is the difference between the two 
metallicity values. The assumption of constant age is sup- 
ported by the fact that we do not observe any age gradient 
in the stellar population (see Fig.jS]). We found that the dis- 
crepancy between observed and predicted g'— i', g'— r', r'— i' 
values is small but it becomes statistically significant for the 
B — I colour index. 

In summary, both galaxies demonstrate the observed 
trend typical of luminous elliptical galax ies, with the galax y 
centre redder than the outer regions (e.g. lFranx et al." 19891 ). 
The SSP predictions for the Sloan colours match the obser- 



NGC 1407 


Profile range 


Obs grad 


SSP grad 




(arcsec) 


(mag dex~^) 


(mag dcx^^) 


B - I 


0.1 - 100 


-0.16 ± 0.01 






1 - 40 


-0.05 ± 0.01 


-0.16 ± 0.01 


g' _ i' 


5 - 100 


-0.12 ± 0.01 






5-40 


-0.10 ± 0.01 


-0.14 ± 0.01 


g' - r' 


5 - 100 


-0.06 ± 0.01 






5-40 


-0.08 ± 0.01 


-0.10 ± 0.01 


r' — i' 


5 - 100 


-0.06 ± 0.01 






5-40 


-0.04 ± 0.01 


-0.04 ± 0.01 


NGC 1400 


Profile range 


Obs grad 


SSP grad 




(arcsec) 


(mag dex~^) 


(mag dex~^) 


g' - i' 


5 - 100 


-0.22 ± 0.01 






5-36 


-0.18 ± 0.01 


-0.15 ± 0.01 


g' - r' 


5 - 100 


-0.11 ± 0.01 






5-36 


-0.12 ± 0.01 


-0.12 ± 0.01 


r' — i' 


5 - 100 


-0.11 ± 0.01 






5-36 


-0.06 ± 0.01 


-0.04 ± 0.01 



Table 5. Observed and predicted colour radial gradients for 
NGC 1407 and NGC 1400. 

vations if ad hoc offsets are applied. On the other hand, we 
found a discrepancy in the slope between the NGC 1407 
observed and predicted radial profile for the B— I colour. 



6 DISCUSSION 

To summarise, we have performed a detailed spectropho- 
tometric study of the early-type galaxies NGC 1407 and 
NGC 1400. From spectral Lick/IDS line-strength indices 
we have measured radial ages, metallicities and abundance 
ratios. Results from high resolution wide-field multi-band 
imaging data (Paper I) have been used to derive radial 
colour profiles. 
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Figure 4. Observed (open circles) and predicted (filled black squares) colour profiles for NGC 1407 (left panel) and NGC 1400 (right 
panel). Colours arc predicted from the measured ages and metallicities using the ?) SSP models. The error bar is the median value for 
the predicted parameters. Offsets have been applied to the predicted data to match our results. 



NGC 1407 


{AC/A[Z/H])obs 


{AC/A[Z/H])ssp 


B - I 


-0.13 ± 0.02 


-0.40 ± 0.01 


g' - i' 


-0.26 ± 0.02 


-0.35 ± 0.01 


g' - r' 


-0.21 ± 0.02 


-0.25 ± 0.01 


r' — i' 


-0.11 ± 0.02 


-0.10 ± 0.01 


NGC 1400 


{AC/A[Z/H])obs 


{AC/A[Z/H])ssp 


g' - i' 


-0.38 ± 0.02 


-0.36 ± 0.01 


g' - r' 


-0.26 ± 0.02 


-0.28 ± 0.01 


r' — i' 


-0.13 ± 0.02 


-0.09 ± 0.01 



Table 6. Observed and predicted 
details. 



AC 
AlZ/H] 



ratio. See Section 5 for 



The lack of radial variation in the luminosity-weighted 
mean age and a-element profiles suggest that the central 
and external galaxy regions formed in a single burst of star 
formation on a timescale of ~ 1 Gyr and then evolved pas- 
sively thereafter. The old central age values, 12.0 ± 1.1 Gyr 
(NGC 1407) and 13.8 ±1.1 Gyr (NGC 1400), constrain the 
redshift of the star forming episode to z ^ 5. Thus neither 
NGC 1407 nor NGC 1400 have experienced any recent star 
formation, at least in a burst strong enough to significantly 
affect the observed integrated light. Donas ct al. (2006) 
combined GALEX UV and optical photometry to study 
the colours of early-type galaxies. The sample included 
NGC 1407 and NGC 1400. UV colours are more sensitive 
to small bursts of star formation. They found the two galax- 
ies to be old with no evidence of residual star formation 
activity. A sim ilar result, also using GALEX, was found by 
lYi et al. (20051 ). The uniformity of the age profiles implies 
that the galaxies have not experienced gaseous mergers since 
formation at redshift z ^ 5. However, we can not preclude 



dissipationless merger events where starbursts are not trig- 
gered. A signature of peculiar dynamics associated with a 
merger event, might be represented by the possible kine- 
matically decoupled core detected in NGC 1407 (Paper I). 
However, the detection is uncertain and possibly due to a 
misalignment of the slit respect the nucleus of the galaxy. 
Moreover, the results of our galaxy isophote modelling show 
no evidence for recent mergers (see Paper I for details). 

The formation of the bulk of the stars at high redshift is 
in concordance with the predictions of the revised monolithic 
collapse models, whereas hierarchical models tend to place 
the formation of elliptical galaxy stars at lower redshift, 
e.g. z < 2. Neverthele ss, the recent numerical simulations of 
iDe Lucia et aL (20061 ) pointed out the distinction between 
formation time of the stars in the progenitors, and the as- 
sembly time of the final ellip tical galaxy, with particula r em- 
phasis on BCGs and BGGS l|De Lucia fc Blaizot 20071 ). Our 
results are in general agreement with their predictions for 
the peak of the star formation at redshift z ~ 5 ('^ 12 Gyr) 
and the formation of ~ 70 per cent of the stars by z ~ 2 
(~ 10 Gyr). However, we find that our results do not repro- 
duce the predicted "bursty" behaviour for the star formation 
history, due to the assumption that bulge stars form during 
merger-induced bursts. We also find a discrepancy in the du- 
ration of the star formation episode (see below), predicted 
to last for several Gyr. 

Recently, iPipino et al. (20071 ) has modelled the [a/Fe] 
radial trend by an outside-in dissipative formation. It is a 
consequence of the interplay between the velocity of the 
Q-enhanced radial flows and the intensity and duration of 
the star formation process at any radius. Thus, the ob- 
served flatness of the NGC 1407 and NGC 1400 [a/Fe] gra- 
dients should be considered a consequence of the faster in- 
flow rate of the a-enhanced gas with respect to the rate 
at which the gas is consumed by star formation, result- 
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ing in a pile-up of gas in the centre. I Thomas et al. (20051 ) 
suggested that the [a/Fe] ratio is proportional to the log- 
arithm of the timescale of a star forming episode, which 
is the FWHM of a Gaus sian star formation ra te. The link 
was previously shown by iThomas et al. and it is in 

ve ry good agreement with the chemodynamical simulations 
of IPipino fc Matteucci f20(M ). They suggest that a ratio 
[a/Fe] = 0.2 for a composite stellar population requires for- 
mation timescales < 1 Gyr. The values that we found suggest 
star formation episodes of a time-duration less than or equal 
to 1 Gyr (star formation rates 100 MQ/year). 

We next turn our attention to the strong metallic- 
ity gradients present in NGC 1407 and NGC 1400. The 
study of such gradients provide an additional diagnostic 
into the mechanism of galaxy formation and evolution, and 
a key constraint to discern between the different theoreti- 
cal models. The measured [Z/H] gradient is —0.38 ± 0.04 
for NGC 1407 and -0.47 ± 0.04 for NGC 1400. Such 
strong gradients are inconsistent with the predictions of 
iKobavashi (2004 ) for any of the dissipative or dissipationless 
merger models. In these models, small or null gradients are 
predicted due to the "washing-out" of the ambient gradients 
in the turbulent mixing of merging events. If star forma- 
tion happens the central metallicity would increase slightly, 
creating a shallow gradient. The gradients in our galaxies 
are steeper than are achieved in either dissipative or dis- 
sipationless merger models. In addition, we can preclude 
that the galaxies have experienced a recent major merger 
due to the lack of tidal tails, shells and plumes expected 
to persist for ~ 3 Gyr after a major merger. At the same 
time, we can exclude the alternative of a minor merger: i.e. 
th e accretion of a s mall, gaseous satellite galaxy (as defined 
bv IKobavashi 20041 ). In this scenario, the satellite galaxy is 
expected to contain low- metallicity gas and therefore high 
central metallicities are hard to achieve. On the other hand, 
steep metallicity gradients are a natural outcome of the re- 
vised dissipative collapse models. The steepness of the slope 
is a consequence of gas dissipation towards inner radii and 
the time-delay in the occurrence of galactic winds due to the 
difference in the galaxy local binding energy. 

ICenarro et al. f2007l ) performed a Lick/IDS spectral 
analysis of a small sample of globular clusters in NGC 1407. 
We calculated the mean value of the derived stel- 
lar populations for the red (metal-rich) globular clus- 
ter subpopulation, which have been shown to be most 
closely associated with the underlyi n g galaxy light (e.g . 
iForbes et al. 20041 : iForbes et al. 20061 : iNorris et al. 20071 ). 
We find: mean age = 10.2 ±2.1 Gyr; [Z/H] = -0.17 ± 0.08; 
[a/Fe] = -1-0.27 ± 0.07. In order to compare these values 
at a similar galactocentric radius, we slightly extrapolate 
our stellar population profiles (Fig. [5} to 1 effective radius 
(the average value for the globular cluster locations) find- 
ing: mean age = 12.5 ± 1.1 Gyr; [Z/H] = -0.21 ± 0.02; 
[a/Fe] = -1-0.45 ± 0.08. Thus, the stellar populations of 
the metal-rich globular clusters and underlying galaxy light, 
at a similar position in the galaxy, are similar to each 
other. This suggests a common formation history for the 
metal-rich globular clu sters and the galaxy spheroid (see 
iBrodie fc Strader 20061 ). 



7 CONCLUSIONS 

We find that the stellar population of NGC 1407 is uniformly 
old, 12.0 ± 1.1 Gyr, with a supersolar degree of a-elements 
enhancement, 0.41 ± 0.05. We also measure a steep negative 
metallicity gradient of —0.38 ± 0.04 and a central metallic- 
ity of 0.29 ±0.08. The stellar population of NGC 1400 is 
found to be constantly old, 13.8 ± 1.1 Gyr, and largely en- 
riched by a-elements, 0.33 ± 0.04. The central metallicity is 
0.25 ± 0.06 and it steeply declines towards outer radii with 
a gradient of -0.47 ± 0.04. 

The results of this work (and Paper I) suggest a similar 
star formation and evolutionary history for the two galaxies. 
The outlined scenario is compatible with the revised version 
of the monolithic collapse model of galaxy formation and 
evolution. The two galaxies formed the bulk of their stars 
at z > 5. The star forming episode lasted for no longer than 
1 Gyr with a star formation rate ^ 100 Mq /year. NGC 1407 
might have experienced an ancient merger, as inferred by 
the possible kinematically decoupled core (Paper I); never- 
theless, the detection is uncertain and potentially originated 
by a misalignment of the slit with respect to the centre of 
the galaxy, as suggested by the ha radial profile. In Paper I 
we have also claimed that NGC 1400 has not interacted with 
NGC 1407 or the group intergalactic medium. The fiatness of 
the a-element radial profiles and the strong negative metal- 
licity gradients are consistent with an outside-in formation 
scenario, where supernova-driven galactic winds played a 
fundamental role and shaped the slope of the gradients. The 
uniform radial profile of old ages, steep metallicity gradients 
and lack of signature of recent mergers indicate that the 
galaxies formed stars quickly some ~ 12 Cry ago and then 
have evolved quiescently without any significant interaction 
ever since. 
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APPENDIX A: CALIBRATIONS TO THE 
LICK/IDS SYSTEM 

The five following Lick/IDS and spectrophotometric stan- 
dard stars were observed: HD4628 (K2), HD4656 (K5), 
HD69267 (K4), HD83618 (K2.5), HD97907 (K3). 

After the correction for recession velocity redshift, the 
standard stars spectra were broadened to the Lick/IDS res- 
olution. The procedure creates a resultant spectrum of total 
broadening: 

o-G = o-v + cri , (Al) 

where ay is the velocity dispersion and a\ the ap- 
propriate Lick/IDS calibration resolution. The Lick/IDS 
spectral r esolution for each index is estimated from 
Fig. 7 of IWorthev fc Ottaviani fl997h . as in Table 2 of 
IProctor fc Sansom f2002t) . The Lick/IDS calibration reso- 
lution: 

= + al , (A2) 
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Index 


Unit 


Lick/IDS Offset 


Error in Mean 




A 


-0.229 


0.120 


v_/ CLrtrr 


A 


-0.054 


0.171 


Fe4531 


A 


-0.342 


0.103 


C4668 


A 


-0.598 


0.187 




A 


-0.041 


0.113 


Fe5015 


A 


-0.357 


0.185 


Mgi 


mag 


0.019 


0.007 


Mg2 


mag 


0.018 


0.006 


Mgt 


A 


0.039 


0.195 


Fe5270 


A 


-0.412 


0.147 


Fe5335 


A 


-0.274 


0.215 


Fo5406 


A 


-0.302 


0.059 


Fe5709 


A 


0.031 


0.065 


Fe5782 


A 


0.026 


0.067 


NaD 


A 


-0.239 


0.138 


TiOi 


mag 


0.006 


0.003 



NGC 1407 



Table Al. Offsets between our index measurements and pub- 
lished values of the five Lick/IDS and spectrophotometric stan- 
dard stars. 



is reached taking into account the instrumental broadening 
of our data erf and convolving the spectrum with a Gaussian 
of width ag. At this stage Lick/IDS indices were measured 
from the standard stars spectra and offsets between our in- 
dex measurements and published values for the Lick/IDS 
standard stars estimated (Table \KT\ . 

The Fourier quotient technique within IRAF software 
package was used to estimate recession velocity and velocity 
dispersion for each aperture spectrum. The kinematic mea- 
surements are in good agreement with those obtained with 
van der Marel's technique (Paper I). We calibrated adopting 
the same technique for the standard stars. The method is 
applicable when the measured velocity dispersion satisfies: 



2 I 2 , 



(A3) 



When the opposite condition is found, the aperture spectra 
are left un-broadened. Lick/IDS indices were then measured 
and index correction factors applied. 

Th index corre ction factors were estima ted using the 
method described in lProctor fc Sansom (20021) (Appendices 
Al and A2). Briefly, the Lick/IDS indices were measured in 
the 5 Lick/IDS standard stars after convolving the spectra 
with a series of Gaussians of known widths. Consequently, 
the appropriate correction factor for a specific index i was 
obtained by a polynomial fit of order 3 



Ci = xo + xiac + X20% -f x-iOc 



where ac = ^J a'j + a 



(A4) 

to the correction factors of the 
index i (average of the values from the 5 standard stars), at 
each of the Gaussian broadenings. 

The offsets evaluated from the Lick/IDS library stars 
were applied to the Lick/IDS indices measured from the ID 
galaxy spectra. Index errors arise from two types of errors 
which are added in quadrature: i.e. random statistical errors 
due by Poisson noise and uncertainties in recession veloc- 
ity and velocity dispersion and systematic errors originated 
from the conversion onto the Lick/IDS system, quantified by 
the error on the mean of the offsets. Systematic and random 
errors are of the same order. 



Index 




A/ndex 
A log r 






Indexr 






Indexo 




-u. 


,65 


± 


0, 


1 A 




,40 


± 


0, 


10 


r. 
0. 


,21 


± 


0.46 




n 

-u. 


.27 


± 


0, 






,49 


± 


0, 


1 n 


^ 


,86 


± 


0.29 


Fe4531 


n 

-u. 


,37 


± 


0, 


1 n 


2 


,77 


± 


0, 


12 


Q 
O. 


,24 


± 


0.35 




_2 


,50 


± 


0, 






,55 


± 


0, 


1 s 


7 


,75 


± 


0.53 


H,3 


-0, 


,01 


± 


0, 


,07 


1, 


,28 


± 


0, 


,08 


1, 


,29 


± 


0.22 


Fe5015 


-1, 


,83 


± 


0, 


14 


3, 


,44 


± 


0, 


16 


5, 


,84 


± 


0.48 


Mgi 


-0, 


,04 


± 


0, 


,00 


0, 


,12 


± 


0, 


,00 


0, 


,18 


± 


0.01 


Mg2 


-0, 


,08 


± 


0, 


,00 


0, 


,24 


± 


0, 


,00 


0, 


,34 


± 


0.01 


Mg6 


-1, 


,08 


± 


0, 


,09 


3, 


,71 


± 


0, 


,10 


5, 


,12 


± 


0.30 


Fe5270 


-0, 


,81 
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0, 


,08 


1, 


,64 
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0, 


,09 


2, 


,68 
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Fe5335 
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1, 


,49 
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NaD 
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,05 
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,07 


2, 


,41 
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,07 
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,01 
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Table Bl. Results of linear regression analysis to NGC 1407 
Lick/IDS indices. Logarithmic radial gradients, values at the ef- 
fective radius (measured by the intercept at r = re) and central 
index values (average of the values extracted in all the apertures 
within a radius of re/8) are given. 



NGC 1400 


Index 


A log r 


IndeXr^ 


Indexo 



Fe4383 


-1.16 ± 0.17 


3, 


,74 


± 


0, 


14 


5, 


,10 


± 


0, 


,45 


Ca4455 


-0.43 ± 0.10 


1, 


,22 


± 


0, 


,08 


1, 


,74 


± 


0, 


,29 


Fe4531 


-0.80 ± 0.13 


2, 


,49 


± 


0, 


11 


3, 


,46 


± 


0, 


,34 


C4668 


-2.82 ± 0.19 


4, 


,19 


± 


0, 


16 


7, 


,47 


± 


0, 


,52 




0.11 ± 0.08 


1, 


,39 
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0, 


,07 


1, 


,26 


± 


0, 


,21 


Fe5015 


-1.65 ± 0.17 


3, 


,79 


± 


0, 


14 


5, 


,63 


± 


0, 


,47 


Mgi 


-0.06 ± 0.00 


0, 


,11 


± 


0, 


,00 


0, 


,17 


± 


0, 


,01 


Mg2 


-0.10 ± 0.00 
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,21 
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0, 


,00 
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,33 
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0, 


,01 


Mgt 


-1.69 ± 0.10 
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,43 
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,08 
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,36 
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,30 


Fe5270 


-0.85 ± 0.09 
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,81 
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,26 


Fe5335 


-1.16 ± 0.12 
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,94 
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,36 


Fe5406 


-0.62 ± 0.08 


1, 


,17 
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0, 


,06 


1, 


,90 


± 


0, 


,22 


Fe5709 


-0.07 ± 0.05 


0, 


,83 


± 


0, 


,05 


0, 


,91 


± 


0, 


,15 


Fe5782 


-0.43 ± 0.05 


0, 


,58 


± 


0, 


,05 


1, 


,01 


± 


0, 


,15 


NaD 


-2.52 ± 0.08 


2, 


,49 


± 


0, 


,06 


5, 


,37 


± 


0, 


,22 


TiOi 


-0.01 ± 0.00 


0, 


,03 


± 


0, 


,00 


0, 


,05 


± 


0, 


,00 



Table B2. Results of linear regression analysis to NGC 1400 
Lick/IDS indices. Logarithmic radial gradients, values at the ef- 
fective radius (measured by the intercept at r = re) and central 
index values (average of the values extracted in all the apertures 
within a radius of re/8) are given. 



APPENDIX B: LICK/IDS LINE-STRENGTH 
RADIAL PROFILES 

Logarithmic gradients of the indices were measured by ap- 
plying a linear regression analysis to the measured index 
values: 



Indexij-) — IndeXr 



AIndex , r 

+ -n log — 

A log r Te 



(Bl) 
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where r is the galactocentric radius, re is the effective ra- 
dius of the galaxy and IndeXr^ is the index value defined at 
the effective radius. We accounted for observational errors 
introducing the weights = cr~^, where o"i are the errors 
associated to measured values, yi, of a specific Lick/IDS in- 
dex at the position Xi (= log(ri/re)). The former equation 
is then equivalent to y = A + Bx, where 

n n n n 

A = ujixf ^ iUtyt - ^ WiXi ^ LOiXiyi)Ar^ (B2) 

i—l i—l i — l i—l 

and 

n n 

B = LUi ^ ujiXtyi 

i=l i=l 

where 

n n n 

A. = '^^uji'^^ujixl - i^^uJiXif . (B4) 

i — l i—l i — l 

The uncertainties on A and B are: 

n 

(^a;OA-i. (B5) 

i = l 

Tables |Bl](NGC 1407) and|B2](NGC 1400) report the 
results of the linear regression analysis. We derived logarith- 
mic radial gradients for all indices and the indices values at 
the effective radius. Central index values were estimated by 
averaging the values extracted in all the apertures within a 
radius of re/8. 



n n 

— '^^ujiXi'^^ujiyi)A~^ (B3) 

i = l 1=1 



\ 



UJiX. 

=1 



\ 
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Figure Bl. Spatially resolved NGC 1407 Lick/IDS line-strength indices radial profiles. Filled red squares represent values from the 
aperture spectra extracted along the observing axis in the negative binning direction; open circles are the same in the positive binning 
direction (see text for details). The approximate extent of the seeing disc is identified by a dashed line. The error bars express the median 
value in the errors data set in the derived line-strength indices. Gradients (black thick lines) are measured excluding the points inside 
the region affected by the seeing. The quantities are measured in logarithmic radius scaled with the effective radius of the galaxy for an 
easy comparison to literature data. 
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Figure B2. Spatially resolved NGC 1400 Lick/IDS line-strength indices radial profiles. Symbols and description are as in Fig. IBll 



